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Abstract The size distributions of 2D and 3D Voronoi cells and of cells 
of V p (2, 3), — 2D cut of 3D Voronoi diagram — are explored, with the single- 
parameter (re-scaled) gamma distribution playing a central role in the an- 
alytical fitting. Observational evidence for a cellular universe is briefly re- 
viewed. A simulated V p (2, 3) map with galaxies lying on the cell boundaries 
is constructed to compare, as regards general appearance, with the observed 
CFA map of galaxies and voids, the parameters of the simulation being so 
chosen as to reproduce the largest observed void size. 

Key words: surveys ; galaxies xlusters: general ; (Cosmology:) large-scale 
structure of Universe 



1 INTRODUCTION 



The applica t ions 



with 



of the Voronoi diagrams (see 



Vorono: 



1 190(3)) m astrophysics started 



Kianel l)l96fil) where the size distribution in ID as given by random seeds was theo- 

also derived , performing a Monte Carlo 



retically deduced in a rigorous way. 



Kiang 



experiment, the area distribution in 2D and volume distribution in 3D. The idea that area 
and volume distributions follow a gamma-v ariate with argument 4 and 6 respectively was 



later reported as "Kiang's conjecture", see 



Okabe et al 



(1992). The application of the 



Icke fc van de Wevgaert 



Voron oi diagrams to the distribution of galaxies started with 
11987), where a seq uential cluste ring process was adopted in order to insert the ini- 
tial seeds. Later on iPierrd l|1990[) introduced a general algorithm for simulating one- 



dimcnsional lines of sight through a cellular universe . The large microwave background 
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temperature anisotropies over angular scales up to one degree were calcula ted u sing a 



Barrow fc Coles! ( 1990 ) and 



Coles 



Voron oi model for large-scale structure formation in 

(1221). 

The possibility to explain the CFA slices using a fractal distribution of seeds and 
insert ing the galaxies on the faces of the irregular polyhedron was explored by 

<iaai) ■ 

A d etailed Monte Carlo s imula tion of pencil beam-like redshift surveys was carried 
out by 



Zaninett 



Sub ba Rao fe Szalavl (|1992j) : they found that the probability of finding regularity 
varies from 3 to 15 percent depending on the details of the models. 



Another Monte Carlo study was carried out by van de Wevgaert fc Bahull (1,994) 
where three different distributions of nuclei were adopted in order to perform extensive 
statistical analysis of several geometrical aspects of three dimensional Voronoi tessella- 
tion. A new way of partitioning the sp ace into cells characterised by the shape of rhombic 

2003) ; the applic ation is done to the CfA 



dodecahedron has been introduced in 



Kiane 



catalogue and to the IRAS/PS Cz catalogue 



proach has been introduced in 



Kiang ct al. (2004). The void hierarchy ap- 



Shcth fe van de Wevga ert ( 2004) and it explains how 



large-scale structures are function of two parameters, one of which reflects the dynamics 
of void formation, and the other the formation of collapsed objects. 

From the point of view of the astronomical observations a few papers that point 
toward the cellular structure of the universe are now briefly reviewed. On analysing the 



Broadhurst et al 



data f rom four distinct surveys at the north and south Galactic poles 
( 1990) found an apparent regularity in the galaxy distribution with a characteristic scale 
of 128 Mpc. The astronomers t hat analysed the maps of the galaxy distribution up t o 



cz=15000 Km/s , see for example 



de Lapparent et al 



( 1988) and 



Geller fe Huchra 



(1989), 



found large coherent structures: the largest void fo und having a diame ter of 5000 km/s. 



Great advances in the observational data , see 



(1996) and 



Shectman et al 



Folkes et al 



(1999) 



Ratcliffe et al 



11996), brought the limits of the observations at cz=60000 
Km/s and confirmed the existence of voids in the distribution of galaxies. 

The distribution of clusters in rich super-clusters is not isotropic: it is periodic along 
a c ubic lattice appro ximately aligned with the super-galactic coordinates, see for exam- 



ple 



Saar et al 



(2002). 



The Voron o i dia gram s are also used to 



El- Ad fc Piranl (|l997(l and 



Ramella et al 



process the astr onomical data , see 



( 2001 ) . As an example 



Ramella et al 



(2001) 



implemented a Voronoi Galaxy Cluster Finder that uses galaxy positions and magnitudes 
to find clusters and determine their main features: size, richness and contrast above the 
background. 

The starting point is to consider a series of explosions that start at the same time in a 
homogeneous space. The shells connected with the explosions meet on a 3D network given 
by the nested irregular polyhedron. From an astrophysical point of view this network 
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can be realized by a set o f primordial explosions, see 



Charlton fc Schramm (1989) and 



Zaninetti fc Ferrard l|1990|) . described by the Sedov solution in the adiabatic phase: 



*(*) = (?— ) 1/5 = 12.49M P c(^) 1/5 , (1) 

M 1/)' V 7l_ 7 ' 

where t represents the time , E is the energy injected in the explosion , p is the density 
of matter , p = nm , n is number of particles per unit volume , m=1.4mH , wh is the 
mass of the hydrogen , t„ = , E 6i = T?F #^ and n_ 7 = 10 -y part [* lcs cm - a ■ 

The already cited works leave the following questions unanswered or partially an- 
swered. 

Is the "Kiang's conjecture" applicable in 2D and 3D environment with physical pa- 
rameters near to those of the galaxies ? 

What is the probability density function that characterises the two-dimensional sec- 
tional area connected with 3D Voronoi diagrams ? 

Can the averaged area connected with the voids in the distribution of galaxies visible 
on the CFA2 slices be guessed from the theory ? 

Can the number of theoretical voids in the distribution of galaxies in a sphere of 
radius equal to that of the CFA2 slices be deduced theoretically ? 

In order to answer these questions the coefficient of the gamma-variate that charac- 
terises the area-distribution of a 2D Voronoi diagram and the volume-distribution in 3D 
were derived in § 12.51 and in § 13 . 31 respectively. 

The coefficient of the gamma-variate that characterises the sectional area of a 3D 
Voronoi diagram as well some characteristics of the voids in the distribution of galaxies 
were derived in § 13.11 

The observed large scale structures of galaxies are classified as CFA slices , LCRS 
slices or pencil beam surveys : they are simulated in § 01 

The number of seeds necessary to produce a theoretical framework comparable to the 
observed one was computed in § 14.1.11 . 



2 THE PRELIMINARIES 

The type of adopted lattice , the importance of setting properly the boundary conditions, 
the type of seeds that generates the polygons/polyhedron , the concept of unitarian area 
and volume and a first two dimensional scan are now introduced. 



2.1 The adopted lattice 

Our method considers a 2D and a 3D lattice made of pixels 2 and pixels 3 points : present in 
this lattice are N s seeds generated according to a random process. All the computations 
are usually performed on this mathematical lattice; the conversion to the physical lattice 
is obtained by multiplying the unit by 5 = pi ^l_ 1 , where side is the length of the 
square/cube expressed in the physical unit adopted. 
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2.2 Boundary conditions 

In order to minimise boundary effects introduced by those polygons/polyhedron that 
cross the square/cubic boundary, we amplify the area/cube in which the seeds are inserted 
by a factor amplify . Therefore the N seeds are inserted in an area/volume that is pixels 2 x 
amplify or pixels 3 x amplify , which is bigger than the box over which we perform the 
scanning; amplify is generally taken to be equal to 1.2 . This procedure inserts periodic 
boundary conditions to our square/cube. The number of seeds that fall in the area/cube 
is N s with N s < N . In order to avoid computing incomplete area/volumes we select the 
cells that do not intersect the square/cubic boundary. This is obtained by selecting the 
cells that belong to seeds which are comprised in an area/volume that is select times 
smaller than pixels 2 or pixels 3 ; of course select is smaller than one and varies between 
0.1 and 0.5. 



2.3 The seeds 



The points are generated independently on the X and Y axis in 2D ( adding the Z axis 
in 3D ) through a subroutine that returns a pseudo-random real number taken from a 
uniform distribution between and 1 . For practical purposes we used the subroutine 

(1992) . 



RANI as described in 



Press et al 



2.4 The adopted units 

In order to deal with quantities of the order of one we divide the obtained area/volume in 
units of pixels 2 / pixels 3 by the expected unitarian area/volume , ua and uy respectively 

ua = pixels 2 x amplify /N 
uy = pixels 3 x amplify /N 

This operation represents a first normalisation. The expected unitarian quantities can 
also be expressed in physical units, u A and Uy 

u p A — side 2 x amplify I N 
Uy = side 3 x amplify /N . 

Special attention should be paid when the area is computed like a cut of a 3D network; 
this case is named V^,(2,3) , see Sec l3.ll . The unitarian area ma(2,3) is expected to be 

UA(2,3) = Uy , (2) 

and the physical counterpart 

< (2 ,3) = «) 2/3 ■ (3) 
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Fig. 1 The Voronoi-diagram in 2D when random seeds are used. The param- 
eters are pixels — 800 , N = 180 ^amplify — 1.2 , side = 2 x 16000 Km/sec , 
select = 0.5 



2.5 The two dimensional scan 

A lattice made of (pixel) 2 points is considered and a typical run using the seeds as given by 
a random process is visualised in F igure ffl Once the histogram of the area is obtained, we 



can fit it , following 



Kiand (|l£66), with the following one parameter probability density 



function , in the following pdf: 

H(x;c) = -f-(ca;) c ~ 1 exp(-cx) , (4) 

r( c ) 

where < x < oo , c >0 and T(c) is the gamma function with argument c. This pdf 
is characterised by /«=1 and a 2 —l/c. The value of c is obtained from the method of the 
matching moments 

1 _ n-l 

° 2 ElU(zi-i) 2 ' [ ) 

The data should be normalised in order to have x =1. The frequency histogram and 
the relative best fit through the gamma- variate are plotted in Figure |21 ■ The captions 
of Figure |2 report also the following quantities expressed in normalised area units (see 
formula iJU)) A, A max , A m i n , respectively the averaged, the maximum and the minimum 
of the area-sample and \ 2 > that represents the goodness of fit. 

The value of c in the case of random seeds is connected with the "Kiang's conjecture" 
that in 2D means c=4 ; the numerical evaluations give similar values , here c = 4.38 is 
found. 
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Fig. 2 Histogram (step-diagram) of area distribution in 2D with a superposition 
of the fitting line (the gamma- variate); input parameters as in Figure^ c=4.38, 
number of bins =10 , x 2 =12.67 , A = 0.98 , A max = 2.19 , A mm = 0.17 . 



3 THE 3D CASE 

In order to make a comparison with the astronomical observations the tessellation in K 3 
is firstly analysed through a planar section and the distribution of volume is numerically 
derived. 



3.1 The 2D cut 

We now work on a 3D lattice Lk. m ,n of pixels 3 elements . Given a section of the cube 
(characterised , for example, by k = p " ei ) the various Vi (the volume belonging to the 
seed i) may or may not cross the little cubes belonging to the two dimensional lattice . 
Following the nomenclature introduced by 



Okabeetal 



1 1992J) we can call the inter- 



section between a plane and the cube previously described as V p (2, 3). A typical result 
of this 2D sectional operation in the x-y plane can be visualised in Figure , the fre- 
quency histogram and the relative best fit through a gamma-variate pdf of the V p (2, 3) 
distribution are reported in Figure 0] together with the derived value of c. Following the 
hypothesis that the galaxies are distributed on the faces of the irregular polyhedra the 
network of Figure |3 represents the spatial coordinates where the galaxies are. The thick 
edges of Figure |3| represent the intersection between the slice and a face ; we remember 
that the area of intersection increases with cz. 

Due to the great importance that the properties of the sample V p {2,2>) could have 
in the astrophysical applications we derived three samples (all crossing the center) of 
Vp(2,3) on the x-y , x-z and y-z plane. This allow us to find the average value of the 
sample properties and their error, see captions of Figure^ 
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Fig. 3 Portion of the Voronoi-diagram V p (2, 3) when random seeds are used; 
cut on the x-y plane . The parameters are pixels = 800 , N — 1900 , side = 2 
x 16000 Km/sec , amplify = 1.2 and select = 0.1 . 




'■ ■ J , I 1 i i L i ■ l- ■ J_i i J i I i i L i i_L i J_. i_ 

D 0.5 1 , 1.5 2 2.5 3 

unitarian area 



Fig. 4 Histogram (step-diagram) of Vp(2, 3) distribution on an x-y plane with a 
superposition of the fitting line (the gamma-variate) ; parameters as in Figure|3 
c=1.99 ± 0.22, number of bin =10 , x 2 =23.09 , A = 0.78 ± 0.06 , A max = 2.12 
± 0.26 , A min = 0.02 ± 0.01 . 



The mathematical theory of the expected values of characteristics of a typical cell in 
1-dimensional sectional Poisson Voronoi diagram, see 



Okabe et al 



gives 



0.68A~ 2/3 



(6) 
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Table 1 The probability to have n-edges in V p (2, 3). 



seeds \ n 


3 


4 


5 


6 


7 


8 


9 


Okabe et al. d992 , l 


0.063 


0.13 


0.2 


0.22 


0.18 


0.11 


0.05 


random 


0.081 


0.16 


0.24 


0.18 


0.10 


0.16 


0.054 




Longitude 



Fig. 5 The Voronoi-diagram V p (2, 3) in the Hammer- Aitof projection at cz — 
7201 Km/sec. 



where A is the intensity of the Poisson process. In our case A is replaced by the random 
points. Our values of A , see captions in Figure^! are near to the values predicted by the 
mathematical theory. In order to obtain the number of edges/crossed faces we report in 
TablcUt he probability to have n edges for random seeds and the theoretical value as given 



Okabe et al 



1 1992). The derived values of c remember the theoretical distribution of 



ID Voronoi segments in which c=2, see 



KiariE (19GG) 



3.2 Projection on the sphere 

Another way to look at the cross sectional area is a spherical cut characterised by a 
constant value of the distance to the center of the box , in this case expressed in cz units, 
see Figure 0. 



3.3 The statistics of the volume 

In every point-lattice l^k.m,n we compute the nearest seed and we increase by one the 
volume of that seed. The frequency histogram and the relative best fit through gamma- 
variate pdf for the volume distribution is reported in FigureEJ The experimental frequen- 
cies are fitted by a gamma-variat e with c = 5 .5 . This value of c should be compared 



with the value of 6 as deduced by 



ch ange in the generato r of random numbers ( 



by 



Kumar et al 



(1992) 



Kiand \l96t ) and success ively refined in 5.5 due to a 



Kiand |l990)) and with 5.78 as deduced 
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Fig. 6 Histogram (step-diagram) of volume distribution with a superposition 
of the fitting line (the gamma-variate) . Parameters as in Figure but pixels = 
400 : c=5.50, NBIN=10 and x 2 =8.05. 



4 THE SPATIAL DISTRIBUTION OF GALAXIES 

The theory of the sectional area derived in Sec. 18.11 can be the framework that explains 
the existence of voids in the spatial distribution of galaxies. The observational material 
that proves the existence of voids in the distribution of galaxies is briefly reviewed and 
then the number of voids in the distribution of galaxies that characterises the CFA2 slices 
is derived. The calibration of the Voronoi diagrams on the void of largest area measurable 
on the CFA slices allows to simulate the maximum in the galaxy's concentration present 
in the various type of slices versus cz and the 128 Mpc regularity. 



4.1 The CFA2 slices 



The second CFA2 redshift Survey , started in 1984, produced slices showing that 
the spatial distribution of galaxies is not random but distributed on filaments that 
represent the 2D projection of 3D bubbles. We recall that a slice comprises all 
the galaxies with magnitude nib < 16.5 in a strip of 6° wide and about 130° 
long. One of such slice (the so called first CFA strip) is visible at the follow- 
in g l address http:// c fa-ww w.harvard .edu~huchra/zcat/| ; more details can be found 



m 



Geller fc Huchra 



I1989J). The already mentioned slice can be down-loaded from 



http: / /cfa-www.harvard.edu~huchra/zcat/n30.dat/| . The first of such sl ices presents 



Huchra 



(2QQ3j)).The great- 



many voids, the bigger one being 4000-^5000 Km/sec across ( 
est possible attention should be paid to the derivation of the void with maximum area 
^mlx > because is the area that allows the determination of the average extension of the 



voids in the distribution of galaxies. The average value of voids diameter can be derived 
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cz (km/s) 



10* 1 JxtO* 



Fig. 7 Polar plot of the little cubes belonging to a slice 130° long and 6° wide. 
Parameters as in Figure [5] 



from the following proportion 



A r , 



Aob 
si™ 



A 



—robs 

A 



(7) 



where the left hand side refers to the maximum and average value of the simulated 
cross-sectional area named V p (2, 3) (their numerical values are visible in the captions of 
Figure|2l and the right hand side refers to the same quantities on the CFA2 slices. A value 



for the average observed diameter : D obs , is easily found from the previous proportion : 



where D, 



obs 



4500 



D 



Ki 



obs 



Q.QD± X = 2700: 



Km 



sec 



27 Mpc 



(8) 



corresponds to the extension of the maximum void visible on 



the CFA2 slices. The half value of D obs can be equated with equation that gives the 
radius of the explosion from primeval galaxies and the following is obtained: 

E6i tg 2 



,J 64 

n_7 



1.47 



(9) 



This relationship regulates the three basic physical parameters involved in the explosions 
of primeval galaxies. 

The results of the simulation can be represented by a slice similar to that observed 
(a strip of 6° wide and about 130° long) , see Figure [7| 

For a more accurate confrontation between simulation and observations the effect 
due to the distribution in luminosity should be introduced . Here conversely a "scaling" 
algorithm is adopted that is now summarised : 

1. The field of velocity of the observed sample is divided in NBIN intervals equally 
spaced. 

2. In each of these NBIN intervals the number of galaxies NGAL(j) , (j identifies the 
selected interval) , is computed. 

3. The field of velocity of the simulated little cubes belonging to the faces is sampled as 
in point (1). 
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Fig. 8 Polar plot of the little cubes (red points) when the "scaling" algorithm 
is applied. NBIN=15 and other parameters as in Figure 0. 



Fig. 9 Polar plot of the real galaxies (green points) belonging to the second 
CFA2 redshift catalogue. 

4. In each interval of the simulated held of velocity NGAL(j) elements are randomly 
selected . 

5. At the end of this process the number of the little cubes belonging to the faces 
equalises the number and the scaling of the observed galaxies. 

A typical polar plot once the "scaling" algorithm is implemented is reported in Figure |H1 
it should be compared with the observations , see Figure 

4-1-1 The density of voids 

The density of seeds expressed in physical units , px, is the inverse of the physical 
averaged volume, pjy — —r, and therefore 




v 




(10) 
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At the same time the sectional area will be characterised by a maximum physical area , 
A$n ax expressed in physical units, 

■^max = C Amax x w ^(2,3) l (H) 

and the following is easily found 

_ ,C_Amax\3/2 / 1 \ 

PN - [—Tp J ■ (12) 

s±max 

We are now ready to compute the number of sources in a sphere of radius R bs = 
16000 Km/ sec , the same radius that characterises the CFA2 slices. The number of 
voids/seeds in the sphere turns out to be 

1.7 l^{C Amax )^ 

JV I Aobs \3/2 ' 

\ max } 

where ^4„ aa; was identified with A°^ ax . Inserting the dimension of the maximum void as 
deduced in g. IPl we obtain A% s ax = 1.59 10 7 '{Km/ sec) 2 . We now have an expression for 
the number of seeds in the sphere that characterises the CFA2 

./V = 268(CAmax) 3 ^ — 827 . (14) 



5 CONCLUSIONS 

The Voronoi diagrams explain and characterise the voids in the spatial distribution of 
galaxies. The characteristics of a typical cell in a two-dimensional section of a 3D Voronoi 
diagram can be compared with those connected with the voids in the distribution of 
galaxies. The following items turn out to be useful to the astronomer once , the 

maximum area connected with a void, is derived from the astronomical slices: 

— The averaged value of the voids in the distribution of galaxies should be 2741 ± 210 
Km/sec across. 

— The pdf of the area of the voids in the distribution of galaxies should be a gamma- 
variate with argument 1.9. 

— The expected averaged value of the sides of the irregular polygons that characterises 
the voids in the distribution of galaxies should 5 . 

Some important key questions are tentatively addressed to the astronomical commu- 
nity: 

— The maximum area connected with a void should be derived with a great accuracy 
in the various slices. 

— The algorithms of describing polygonal voids from the astronomical observations 
should be developed in order to test the suggested averaged number of sides, 5 , 
as predicted from the Voronoi diagrams. 
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— The pdf of the area connected with the voids should be tentatively computed in order 
to test the goodness of the V p (2,3) predictions. 

The Voronoi diagrams allow also to reformulate the theory of the primordial explo- 
sions because 

— The average diameter of voids between galaxies is function of three parameters : time 
, density and energy , see formula El 

— The galaxies are originated where the primordial shells meet , the faces of the irregular 
polyhedra. 

— The correlation length for galaxies can be identified with the face's thickness that is 
approximately | the radius of the expanding shell. 

Acknowledgements I thank the Smithsonian Astrophysical Observatory and John Huchra for 
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